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Abstract. The purpose of this study was to determine the growth variability of four provenances of Picea abies on experimental 
plots in the Wyszkow Forest District, central Poland. The experiment was established as a system of random blocks with four 
repetitions per block. We selected 48 trees from each provenance and the increment cores were colected from sample trees. 
Standard measurements of the width of annual increments were performed using the WinDendro software. Raw data was then 
indexed and subject to dendroclimatic analyses based on the average monthly temperatures and precipitation of the period from 
1969 to 2012. Furthermore, the COFECHA software was used to check the consistency of the data and to determine the pointer 
years. High data consistency as well as growth variability of particular provenances in response to climatic conditions was 
observed. The results obtained here will allow for an improved selection of populations best suited for growing in the climate of 
central Poland. 
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1. Introduction 

Spruce is one of the main forest-forming species in Poland. 
According to Srodon (1976), this species occurs in two di- 
stinct areas in Poland. Central and northwestern Poland lacks 
natural sites of this species, although it has been artificially in- 
troduced there. This is called the ‘spruceless’ area. The origin 
of the spruceless area has been the subject of many years of 
ongoing discussion, resulting in several theories: 

- It is a natural border separating areas of lowland and 
mountain spruce. 

- It is the southern limit of the range of spruce from the 
northwest, designating the species’ unfinished migration to 
the south (Szafer 1931, Schmidt- Vogt 1977). 

- It was created as the result of deliberate forest manage- 
ment in areas unfavourable for the growth of spruce (Hensel 
1980, Broda 1998). 



The last theory is supported by recent studies on the effects 
of climate on the radial growth of spruce in different prove- 
nances (Koprowski and Zielski 2006; Koprowski 2013). 

One of the aims of numerous provenance studies on the 
adaptive possibilities of spruce, originating from its natural 
range as well as outside of it, was to find a population able 
to grow in the habitat of the spruceless area (Giertych 1972, 
2000; Krupski and Giertych 1996). However, mainly gro- 
wth and qualitative characteristics were analysed. To date, 
few studies on the provenance of spruce were devoted to 
the relation between genetic determinants of radial growth 
and climate (Burczyk and Giertych 1988, 1991; Koprow- 
ski 2008). Based on the results of similar studies on other 
conifer species, it can be assumed that such diversity exists 
(Oleksyn and Fritts 1991; Savva et al. 2002; Eilmann et al. 
2013; Bijak, 2013). Spiecker (2002) indicates a genetic basis 
for growth variation in response to climatic conditions. 
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According to Koprowski (2013), projected climate chan- 
ges will most likely equally affect the growth of spruce 
populations from its natural range, as well as from the spru- 
celess area. The most recent dendroclimatic studies indicate 
that the reaction is on a supra-species (affecting all conifero- 
us species) (Bijak 2013) or supra-regional level (Koprowski 
and Vitas 2010). 

The aim of this study is to verify the hypothesis of the 
existence of differences in dendroclimatic reactions among 
northern Norway spruce provenances growing in an experi- 
mental plot located in central Poland. 

2. Materials and method 

The study material consisted of trees from the Norway 
spruce provenance area of the Wyszkow Forest District, Da- 
lekie Forest Range, division 99i (52°41’24”N 21°33’58”E), 
founded in 1963. Four of eleven spruce provenances from 
northern Poland were chosen; two characterised by the best 
productivity - Kartuzy (provenance 11) and Mestwinowo 
(provenance 8), and two by the worst productivity - Augu- 
stow (provenance 4) and Gorowo llawieckie (provenance 7) 
(Matras et al. 2006; Kowalczyk, unpublished data). These 
provenances also represent populations within the natural 



range of spruce - Augustow and Gorowo, as well as in the 
spruceless area - Kartuzy and Mestwinowo (Fig. 1). 

Pith to bark increment cores at 1.3 m height were taken 
from 192 trees, 48 from each selected provenance. The core 
samples were prepared according to standard techniques and 
then the width of annual rings was measured and individu- 
al sequences were established for each tree with WinDEN- 
DRO. Tree-ring sequences were verified: the accuracy of the 
dating and the homogeneity of the series were confirmed 
using the COFECHA program with the DPL package (Holmes 
1999; Grissino-Mayer 2001). Series with the lowest linear 
correlation coefficients and series incompatible with others 
were discarded from further analysis (a total of 10 series were 
rejected). Then, the chronologies of each provenance were es- 
tablished - the actual and residual - with the Arstan DPL pro- 
gram package (Cook 1985; Cook et al. 1990). A 64-year filter 
was applied in the standardisation process. The similarity be- 
tween chronologies was analysed using convergence coeffi- 
cient GL, t-value and correlation coefficient (a = 0.05). The 
amalgamation method of cluster analysis (Tryon 1939) was 
used to compare the growth patterns of the analysed Norway 
spruce provenances. In constructing hierarchical graphs, the 
Euclidean distance was adopted to calculate the proximity of 
the groups. The distances between clusters were determined 




Figure 1. Location of the analysed provenances and diagram of the experimental plot 
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by the complete linkage method, in which the distance be- 
tween clusters is determined by the greatest of the distances 
between any two objects belonging to different clusters. Cli- 
mate-tree growth relationships were examined for the period 
of 1980-2012 using a mathematical model of the response 
function (RESPO program) (Briffa and Cook 1990). Residual 
chronology parameters of the specific provenances were used 
as growth parameters. Climatic parameters (mean monthly air 
temperature and monthly precipitation totals) were obtained 
from the Meteorological and Hydrological Institute station in 
Warsaw. The range of climate variables included 17 months 
of data, from May of the year preceding the formation of tree 
rings to September, when the ring had completed its growth. 
Pointer years were determined using the WEISER program 



(Gonzalez 2001) for each provenance, as well as for all of 
them together. A compliance threshold of 80% with a mini- 
mum number of 10 trees per provenance was adopted. 

3. Results 

The analyses exhibited a high degree of agreement for 
all growth curves, both for individual sequences represen- 
ting specific provenances, as well as for the chronology of 
the provenances. The degree of similarity of the chronology 
constructed for each provenance was high (Figs 2 and 3; Ta- 
bles 1 and 2). The multiple range test used to compare the 
analysed provenances did not confirm significant differences 
among the provenances. 




Years 



• provenance 4 provenance 7 — — provenance 8 provenance 1 1 




Figure 2. Actual (upper graph) and 
residual (bottom graph) chronologies 
of the analysed spruce provenances. 
Forest districts: Kartuzy (provenance 
11), Mestwinowo (provenance 8), 
Augustow (provenance 4) and Gorowo 
Ilawieckie (provenance 7) 
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Kartuzy Gorowo Augustow Mestwinowo Figure 3. The similarity of growth patterns for spruce provenances based 

on cluster analysis 
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Table 1. Characteristics of the chronologies representing each provenance 



Chronology 



No. of 

Provenance 

trees type 



average 



standard 

deviation 



coefficient of 
variation 



average 

sensitivity 



autocorrelation 



Augustow 




raw data chronology 


2.1 


0.8353 


40% 


0.1718 


0.8009 


49 














residual 


1 


0.2226 


22% 


0.2425 


0.0094 


Gorowo 


48 


raw data chronology 


2.01 


0.8337 


41% 


0.1701 


0.8486 






residual 


1 


0.2437 


24% 


0.2626 


0.1127 


Mestwinowo 




raw data chronology 


2.15 


0.8594 


40% 


0.1836 


0.815 


48 
















residual 


1 


0.2104 


21% 


0.2531 


-0.025 






raw data chronology 


2.12 


0.8651 


41% 


0.1578 


0.8935 


Kartuzy 


47 














residual 


1 


0.2788 


28% 


0.2926 


0.09 



Table 2. The degree of similarity of the chronologies representing each provenance 



GL Coefficient [%]) 





Augustow 


Gorowo 


Mestwinowo 


Kartuzy 


Augustow 


X 


79 


79 


85 


Gorowo 


15.80 


X 


85 


81 


Mestwinowo 


17.34 


17.26 


X 


79 


Kartuzy 


20.01 


14.82 


17.20 


X 



The growth reaction of trees from various provenances to 
climatic conditions was similar. The analysis of the response 
function indicated that the relation between growth and cli- 
mate is similar for each series representing a particular pro- 
venance. Thermal conditions affect radial growth mainly in 
the period preceding the formation of the ring - May, June, 
September and November of the previous year (Fig. 4). Pre- 
cipitation conditions of the growing season, both in the year 
preceding ring formation as well as in the year of growth, 
affect its size (Fig. 5). The dependencies presented of ra- 
dial growth to thermal and precipitation conditions were the 
same for each of the four provenance studies. 

Ten pointer years were designated: five positive ones, that 
is, those in which growth increment was especially high (1995, 
2001, 2004, 2007, 2009), and five negative years, in which 
growth was small (1979, 1983, 1986, 2000, 2005). These were 
common to all provenances. Negative years were characterised 
by low rainfall during the growing season. Climatic conditions 
in positive years did not differ from the long-term data pu- 
blished by the Central Statistical Office (Domanska et al. 2009). 



4. Discussion 

A comparison of the growth patterns of the four Norway 
spruce provenances showed no statistically significant dif- 
ferences regardless of the test used (t-value, sign test-G). 
These results are consistent with those obtained by Burczyk 
and Giertych (1988) and Koprowski (2008). Burczyk and 
Giertych (1991) found that the lack of differences between 
the provenances was due to the methodology of choosing 
sample trees, which gave preference to trees of similar dbh 
values. Note, however, that the present study selected trees 
from the I, II and III Kraft biosocial classes. Burczyk and 
Giertych (1991) also indicate a significant influence of mi- 
crohabitat conditions on different dendroclimatic reactions 
among provenances and on the lack of an influence of mi- 
crohabitat conditions on the relationship between radial in- 
crement and climate. Researchers are divided on the issue 
of whether a reaction in radial growth to extreme climatic 
conditions is genetically determined. Although Burczyk and 
Giertych (1991) found no such relationship for the Norway 
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Figure 4. Coefficients of correlation (bars) and regression (points) of average growth values of spruce provenances and air temperature; 
statistically significant values (a = 0.05) are indicated with black bullets and dark grey bars, py - preceding year 





Figure 5. Coefficients of correlation (bars) and regression (points) of average growth values of spruce provenances and precipitation; 
statistically significant values (a = 0.05) are indicated by grey bullets and black bars, py - preceding year 



spruce, Elmann et al. (2013) did in their investigation of Do- 
uglas fir provenances, which depended on the origin of the 
stock. In their research, provenance trees characterised by 
wide annual increments clearly responded to drought. Thus 
the cited authors associated a higher resistance to drought 
stress with low productivity of a population, and therefore 
narrower increments. Savva et al. (2002) reached interesting 
conclusions in research on Russian Scots pine provenances. 
They found that as the distance of the place of origin of a 
provenance increased from the experimental plot, the im- 
pact of local weather conditions decreased. In our study, the 
distance of specific provenances from the test site was simi- 
lar. For this reason, we anticipated a similar effect of clima- 
tic conditions on the radial increment of all four populations. 

The analysis of the similarities of growth patterns among 
the provenances showed a clear difference between the popu- 
lation from Mestwinowo and the others. However, it should 
be recalled that there were no significant differences between 
the growth patterns of the studied provenances. The research 
of Giertych (2000) allows us to explain this phenomenon by 



the forest management and planting conducted by the Prus- 
sian state in this region at the beginning of the 20th century. 
Populations that are now in Poland were managed in a similar 
manner at that time, regardless of whether they were located 
in the natural range of spruce occurrence or in the spruceless 
area. Thus there is no reason to believe that the direction of 
introduction differed for various populations. This may also 
explain the lack of significant differences in growth patterns in 
this study. Koprowski (2013) found that regardless of whether 
spruce populations were from the natural range or not, the in- 
crement response to climatic conditions was the same. 

The results of the relationship of radial increment to cli- 
matic factors were similar for each analysed provenance. 
Previous studies found no correlation between the radial 
increment of spruce and mean temperature in lune of the 
year preceding the analysed growing season. The negative 
correlation of this characteristic with the temperature in 
July, September and November of the preceding year, while 
not statistically significant, was confirmed in the research 
of other authors (Zielski and Koprowski 2001; Bijak 2009; 
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Rybnicek et al. 2012a, b). The nonspecific radial growth re- 
sponse of spruce to the mean temperature in June may be re- 
lated to climatic conditions, which, according to Koprowski 
and Zielski (2003), affect the growth pattern of this species 
and determine its spatial distribution. 

The analysis of correlation between the radial increment 
of spruce and precipitation confirmed a positive and signifi- 
cant relationship with the rainfall in July and September of 
the preceding year, as well as in June and July of the growth 
year. These results correspond with those of areas with similar 
geographical conditions (Zielski and Koprowski 2001; Ko- 
prowski and Zielski 2006; Bijak 2009) as well as from sites 
located both south (Rybnicek et al. 2012a, b) and north (Vitas 
2004) of the study area. The similar nature of radial growth 
response to climate in different locations casts doubt on the 
hypothesis of the varying influence of environmental factors 
on radial growth in different regions (Makinen et al. 2003). 

In analysing the pointer years selected for the present study 
and of other authors, we obtained the same results for only 
one year - 1979 (Koprowski and Zielski 2002; Vitas 2004). In 
that year, weather conditions negatively influenced the radial 
growth of spruce from the analysed provenances at Wyszkow, 
as well as spruce growing in the Olsztyn Lake District and 
western Lithuania. Both positive and negative characteristics 
of the remaining pointer years were unique to our area. Accor- 
ding to Vitas (2001), differences in the designation of pointer 
years may result from the impact of factors other than climate 
and may be genetically determined. The reasons for the dif- 
ferences in pointer years could also be due to the differences 
in the ages of the analysed stands. From 1980 to 2000, the 
stands analysed in this study were between the ages of 20 and 
40 years, while those analysed by other authors were between 
the ages of 180 and 200 years (Koprowski and Zielski 2002) 
and 100 and 120 years (Vitas 2004). In addition, these authors 
conducted their studies in spruce stands, not experimental 
plots, where due to the small size of an individual plot the 
impact of climatic conditions may differ. The growth rates of 
spruce in some pointer years (1995, 2000, 2001, 2004) of the 
analysed 20-year period cannot be explained by climatic ano- 
malies. Therefore, they could have been related to other stress 
factors reducing radial growth, which occurred only locally. 

5. Conclusions 

- There is no evidence for differences in the dendrocli- 
matic reaction of northern Norway spruce provenances to 
growth at the Wyszkow experimental station. 

- Precipitation uniformly influences the radial growth of 
Norway spruce, irrespective of the continentality gradient, 
while thermal conditions result in a nonspecific growth re- 
sponse related to the location of the stand. 



- Increment in negative pointer years not confirmed by 
other published studies may be the result of factors other 
than climate, such as stand age and growth conditions, 
among others, at the study site. 
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